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The NAVSTAR Global Positioning System (GPS) has provided a quantum leap in 
real-time autonomous navigation capabilities. NASA’s Space Shuttle will be receiving an 
integrated GPS capability in the near future, and the orbiter Endeavour has been equipped 
with a stand-alone GPS receiver. Although much data is available regarding spacecraft 
GPS receiver performance at higher altitudes, little information is available for spacecraft 
at Shuttle altitudes of approximately 400 km where drag and gravity effects are more 
pronounced. 

GPS receiver navigation solution data from Shuttle mission STS-69 was made 
available by NASA and provided an opportunity for evaluating GPS performance in low 
Earth orbit. NASA ground tracking network and Tracking and Data Relay Satellite 
(TDRS) data for this mission provided a reference for comparison. Analysis of the data 
was accomplished using Satellite Tool Kit (STK) for visualization and Matlab routines for 
data comparison. 

GPS navigation solutions were available for approximately 65 percent of the STS- 
69 mission, and they generally coincided with the reference track. Differences between the 
GPS navigation solution state vectors obtained using the Standard Positioning Service 
(SPS) and the reference state vectors produced RMS position differences between the 
data sets of about 1500 m. One sigma position accuracies of 54 m in the vertical direction 
and approximately 1400 m in the downtrack direction were experienced. Velocity vector 
magnitude differences during this period were generally ± lm/s, with a RMS velocity 
difference of less than 9 m/s. One sigma velocity accuracies of approximately 4.2 m/s in 
the vertical direction, 2.3 m/s in the downtrack direction and 1.5 m/s in the crosstrack 
direction were experienced. A firm conclusion regarding Shuttle GPS accuracies could 
not be drawn because all sources of error were not identified. Based on these results GPS 
appears to be an excellent navigation source for Shuttle state vector information; however, 
another navigation source such as INS must be present to provide a check against 
spurious data points and periods of outage. 
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INTRODUCTION 



The NAVSTAR Global Positioning System (GPS) has provided a quantum leap in 
real-time autonomous navigation capabilities. GPS technology has been applied in 
numerous fields and has now become a method for autonomous spacecraft navigation. 
Flight experience with the Topex/Poseidon spacecraft has indicated that very precise 
spacecraft orbit information can be obtained based on GPS measurements. NASA’s 
Space Shuttle will be receiving an integrated GPS capability in the near future, and the 
shuttle Endeavour has been equipped with an operational stand-alone system. Although 
much data is available regarding spacecraft GPS receiver performance at higher altitudes, 
little information is available for spacecraft at Shuttle altitudes of approximately 400 km 
where drag and gravity effects are more pronounced. 

GPS receiver navigation solution data from Shuttle mission STS-69 has been made 
available by NASA and provides an opportunity for evaluating GPS performance in low 
earth orbit. NASA ground tracking network and Tracking and Data Relay Satellite 
(TDRS) data for this mission was also available for comparison. The objective of this 
thesis was to compare Shuttle GPS receiver navigation solutions with the NASA reference 
track. Data comparison routines written in Matlab were used to analyze both sets of state 
vectors. The visualization features of the orbit analysis software package Satellite Tool 
Kit (STK) aided the analysis of the data sets. 

A. GPS BASICS 



1. System Operation 

GPS was developed by the Department of Defense (DoD) as a worldwide satellite- 
based radionavigation system which provides users with position, velocity and time 
information. Using the concept of ranging, the system measures the distance to several 
satellites to determine position. Ranging is dependent on time synchronicity among the 
receiver clock and the satellites. The satellites transmit radio pulses in the form of distinct 
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pseudo-random sequences at specific known times. The receiving equipment recognizes 
these sequences and measures the precise time at which the pulses arrive. Based on the 
transmission and reception times, the time it took the pulses to travel from the satellite to 
the receiver can be determined. Since the speed of light is constant and the travel time of 
the pulses is now known, the range to the satellite can be calculated. 

Determining the range to a particular satellite places the user on a sphere with a 
radius equal to this range. The intersection of two spheres of position is a circle, and the 
intersection of three spheres is two points. Position can be determined using three satellite 
ranges because the ambiguity of the second point can usually be resolved as an 
unreasonable solution. The intersection of four spheres of position produces a single 
point. (Herring, 1996, p. 46) 

Ranging requires that the receiver’s clock be synchronized with the satellite’s 
clock; however, in order to support receivers with less than perfect chronometers, it is 
necessary to apply correction and this requires a fourth satellite as shown in Figure 1.1. 
First, ranges to the four satellites are calculated using the original receiver clock time and 
are called pseudo-ranges. The resultant four spheres of position would intersect in a 
single point if there were no clock error. The receiver clock error is then determined 
algebraically. Because there is only one value of clock error for which the four spheres 
can be made to intersect at a single point, the user’s position solution is thus provided as 
shown in Figure 1.2. (Herring, 1996, pp. 46-47) 

Arranged in six orbital planes with four satellites in each plane, the constellation 
of GPS satellites consists of twenty-four satellites. This ensures that a minimum of four 
satellites are always observable by a user anywhere in the world. The satellites are in 
circular orbits inclined at fifty-five degrees with altitudes of 20,200 km and complete an 
orbit every twelve hours. 



2 




Figure 1.1. Clock error in the receiver causes the range to the GPS satellite to be 
incorrect. The resulting spheres of position (thick lines) do not intersect at a single point. 
By adjusting the receiver clock the ranges are corrected and meet at one point (thin lines). 
In this figure this appears to require three satellites; however, in three dimensions four 
satellites are required. From Herring, 1996, p. 46. 
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2 . 



GPS Accuracy 



GPS provides two levels of service, the Standard Positioning Service (SPS) and 
the Precise Positioning Service (PPS). Utilizing spread spectrum techniques, GPS 
satellites transmit on two L-band frequencies, LI at 1575.42 MHz and L2 at 1227.6 MHz. 
Two spreading functions known as C/A-code, or Coarse/ Acquisition code, and P-code, or 
Precise code, are employed. Currently only the C/A-code signal is available to all GPS 
users and provides the SPS, while the P-code is available to DoD-designated users only 
and is necessary for the PPS. Although user position can be determined using both codes, 
only the P-code is transmitted on both frequencies. This enables users with P-code access 
to perform a dual frequency comparison to compensate for ionospheric effects while C/A 
code users must rely on a less accurate model of the ionosphere. Satellite ephemeris 
(orbital element) data, almanac data for finding other satellites in the constellation, 
atmospheric propagation correction data, satellite clock-bias information, system time and 
satellite status information are provided by the navigation message which is transmitted on 
both frequencies. (U.S. Naval Observatory, 1996, pp. 1-2) 

Full accuracy is denied to non-DoD designated users through the use of Selective 
Availability (SA), which modifies the navigation message. SA reduces accuracy by 
altering the satellite ephemeris data and clock frequency, a procedure which is known as 
dithering. Anti-spoofing protects against false satellite transmissions by encrypting the P- 
code, forming the Y-code. (U.S. Naval Observatory, 1996, p. 2) 

The SPS has a 95 percent probability of providing horizontal positioning accuracy 
within 100 meters, vertical positioning accuracy within 156 meters and timing accuracy 
within 340 nanoseconds. (U. S. Naval Observatory, 1996, p. 1) Typically, one sigma 
positioning accuracy is on the order of 50 meters horizontally and 75 meters vertically 
with a velocity accuracy of 50 cm/s. The specification for PPS requires positioning 
accuracy of 16 meters Spherical Error Probable (SEP) in contrast to the SPS specification 
of 76 meters SEP and a velocity accuracy of 10 cm/s. A positioning accuracy of 10 
meters SEP is usually experienced. (Clynch, 1996) Sources of GPS inaccuracy and their 
impact are shown in Table 1.1. 
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Component 


Size 


Receiver Clock 


1000 km 


Selective Availability 


30 m 


Ionosphere 


1 -30 m 


Atmosphere 


1 - 6 m 


Orbit Error 


1 - 3 m 


Satellite Clock 


1 - 3 m 


Multipath - C/A Code 


0.5 - 150 m 


Multipath - P Code 


0.1-15 m 


Receiver Noise 


0.1 - 1 m 



Table 1.1. GPS Sources of Inaccuracy. FromClynch, 1996. 

GPS system time is maintained by the Composite Clock which incorporates all 
operational Monitor Station and satellite frequency standards. This system time is 
referenced to the U.S. Naval Observatory’s Master Clock within a standard tolerance of 
one microsecond Over the last several years GPS system time has maintained a tolerance 
within a few hundred nanoseconds of the Master Clock. (U.S. Naval Observatory, 1996, 
P-3) 

B. SHUTTLE GPS 

An integrated GPS navigation capability has been developed for the Shuttle. The 
scheduled replacement of TACAN air navigation stations has been the primary motivating 
factor in developing this capability. Currently TACAN provides the primary Shuttle entry 
navigation aid from post blackout to final approach. During final approach through 
landing the Microwave Scanning Beam Landing System (MSBLS) is the primary 
navigation aid. Presently, only the orbiter Endeavour is equipped with GPS equipment. 
The other orbiters will be equipped as they are upgraded. 

Two options were considered for integrating GPS into the Shuttle navigation 
system. The first option was to augment Shuttle entry navigation filters to process GPS 
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measurements. The second option was for the Shuttle navigation system to process state 
vector outputs from the GPS receiver. Either option can be extended to other Shuttle 
mission phases such as ascent and on-orbit operations. The second option of processing 
GPS state vector outputs was selected because it required fewer software changes and had 
the minimum impact on current system operation. Additionally, it facilitated the design 
objective of a phased- development approach which would permit evaluation of GPS 
receiver performance during flight before it was used as the primary entry navigation aid. 
(Kachmar et. al., 1993, pp. 313-317) With respect to developmental testing for STS-69, 
the GPS receiver navigation solution was not available to the orbiter in real time even as a 
back up. Beginning with STS-79, the navigation data will be processed by the orbiter 
computers and real-time state vectors will be available during launch and landing. 

Using the second approach, a GPS receiver state vector can be used to reset the 
baseline navigation state. This navigation state reset procedure is similar to a ground state 
vector uplink with the important difference that the GPS state vector quality assessment is 
done onboard the Shuttle. In order to be selected for guidance, the GPS state vector must 
meet a predefined transient tolerance level. Transients can be induced by the Shuttle 
mission phase environment, the GPS receiver switching between satellites, and periods of 
satellite reacquisition. The GPS receiver can be provided aiding data by the Shuttle’s 
Inertial Measurement Units (IMU) to assist the receiver in determining its position for 
acquisition of the GPS satellites. (Kachmar et. al., 1993, p. 317) 

The Shuttle was designated as a DoD-authorized user of PPS due to the need to 
guarantee a minimum level of entry navigation performance during times of national 
emergency. It was also desired that current off-the-shelf receiver designs be used that 
complied with the GPS Joint Program Office Requirements Document. (Kachmar et. al., 
1993, p. 317) The Collins 3M receiver was selected for installation onboard Endeavour. 
The five channel 3M receiver tracks four satellites simultaneously while reading the 
navigation message for the next satellite to be used by the receiver. It is a P-code 
receiver with position accuracy of 15 meters SEP, velocity accuracy of 10 cm/s, and 
timing accuracy of 40 nanoseconds. (Collins, 1994) (Nuss, e-mail, 1996) Two antennas 
are employed to provide maximum coverage for the receiver. The upper-hemi and lower- 
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hemi antennas are located on the top and bottom of the orbiter respectively as shown in 
Figure 1.3. 



UPPER S-BAND 
HEM) /GPS 




Figure 1.3. Shuttle GPS Antenna Placement. From Carpenter and Hinkel, 1995, p. 6. 



C. STS-69 

Shuttle mission STS-69 lifted off from Pad 39-A at the Kennedy Space Center 
(KSC) on September 7, 1995, at 11:09:00.052 a.m. EDT. Endeavour’s crew of five 
included: Commander David M. Walker, Pilot Kenneth D. Cockrell, Payload Commander 
James S. Voss, Mission Specialist James H. Newman, and Mission Specialist Michael L. 
Gemhardt. The mission featured the second flight of the Wake Shield Facility (WSF), a 
four-meter-diameter saucer-shaped free-flying satellite designed to grow semiconductor 
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films in the ultra-vacuum created in the wake of the satellite as it moves through space. 
The Spartan 201 free-flying astronomy satellite was also deployed and retrieved. 
Additionally, the crew tested assembly techniques for the international space station and 
tested thermal improvements made to space suits used during space walks. On September 
18, 1995 at 7:37:56 am EDT, after 10 days, 20 hours, 28 minutes and 55 seconds, STS-69 
landed at KSC. After a successful mission and traveling over four and a half million miles. 
Endeavour touched down on Runway 33 at the KSC Shuttle Landing Facility. 

The WSF was employed to conduct several GPS experiments sponsored by the 
Texas Space Grant Consortium and was equipped with a TurboRogue GPS receiver 
furnished by the University Corporation for Atmospheric Research (UCAR). University 
of Texas at Austin researchers and Johnson Space Center personnel conducted the 
experiments which were aimed at: assessing relative positioning using the TurboRogue 
receiver on the WSF and Endeavour’s Collins 3M receiver, evaluating precision orbit 
determination in a low-altitude environment (400 km) where perturbations due to 
atmospheric drag and the Earth’s gravity field are more pronounced than for higher 
altitude satellites such as TOPEX/POSEIDON, and determining atmospheric temperature 
profiles using GPS signals passing through the atmosphere. (Schutz et. al., 1995, pp 1-2) 
During STS-69 Endeavour’s 3M receiver operated in a SPS mode. 
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O. REFERENCE FRAMES AND SHUTTLE TRACKING 



A. REFERENCE FRAMES 

Several reference frames were used when comparing GPS navigation solutions 
with the reference track. Shuttle navigation and tracking sources utilize the Aries Mean- 
of-1950 Earth-Centered Inertial reference frame, Greenwich Mean Time and English units. 
GPS navigation solutions provide position data in the WGS 84 Earth-Centered Earth- 
Fixed reference frame using Coordinated Universal Time and metric units. GPS velocity 
data is in a Shuttle Up-East-North frame and in metric untis. 

1. Inertial Reference Frames 

The Aries Mean-of-1950 (M50) and J2000 Earth-Centered Inertial (ECI) reference 
frames are not actually fixed with respect to the mean position of the stars in the vicinity 
of the Sun, or inertial space. These inertial frames are defined by the direction of the 
Earth’s axis of rotation, or celestial pole, and the direction from the Earth to the Sun on 
the first day of spring which is known as the vernal equinox or first point of Aries. The 
celestial pole is affected by gravitational forces, primarily those exerted by the Sun and 
Moon. These forces produce a quasi-conical motion of the mean celestial pole around the 
pole of the ecliptic, or axis of rotation of the Earth’s orbit around the Sun. This lunisolar 
precessional motion has a period of 26,000 years. The motion of the actual celestial pole 
around the mean pole is known as nutation and has a period of 18.6 years. (Larson and 
Wertz, 1992, pp. 94-95) (Seidelmann, 1992, p. 12) The International Astronomical Union 
(IAU) has adopted conventions for calculating general precession (IAU 1976) 
(Seidelmann, 1992, p. 103) and general nutation (IAU 1980) (Seidelmann, 1992, p. Ill); 
however, they do not account for transient pole wander. 

Motion of the vernal equinox subsequently accompanies the motion of the celestial 
pole. As a result, in order to define an ECI reference frame the date which specifies the 
position of the vernal equinox must be established. (Seidelmann, 1992, p. 12) The J2000 
frame is a right-handed Cartesian coordinate system with its origin at the center of the 
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Earth and its X and Z axes pointing towards the mean vernal equinox and mean rotation 
axis of January 1, 2000. The M50 coordinate system uses the mean vernal equinox and 
mean rotation axis of 1950 as its reference directions. A constant transformation matrix 
has been developed to perform the conversion from M50 to J2000 coordinates. 
(Seidelmann, 1992, pp. 184-187) 



2. WGS 84 Reference Frame 

The World Geodetic System 1984 (WGS 84) reference frame is a Conventional 
Terrestrial System (CTS) whose origin is the center of mass of the Earth. The system’s Z- 
axis is parallel to the direction of the Conventional Terrestrial Pole (CTP) defined by the 
Bureau International de l’Heure (BEH). The X-axis is the intersection of the WGS 84 
reference meridian plane, which is parallel to the Zero Meridian defined by the BEH, and 
the CTP’s equatorial plane. The Y-axis completes a right-handed, earth-fixed orthogonal 
coordinate system. This reference frame rotates with the Earth and assumes that the earth 
is rotating at a constant rate around a mean celestial pole, the CTP. Variations from this 
approximation cause the WGS 84 CTS to differ from the true Instantaneous Terrestrial 
System (ITS) which is rotating around an instantaneous pole, the Celestial Ephemeris Pole 
(CEP). The WGS 84 CTS can be related to the J2000 Conventional Inertial System (CIS) 
by Equation 2.1. 



CTS = [A] [B] [C] [D] CIS (2.1) 

The rotation matrices account for: [A] polar motion, [B] sidereal time, [C] 

astronomic nutation, and [D] precession. Matrices B, C and D establish the relationship 
between the CIS and the ITS. The ITS is related to the CTS by matrix A which provides 
the relationship between the CEP and the CTP. (DMA, 1987, pp. 2-1 - 2-3) 
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3. 



Transformation to J2000 



The Shuttle navigation system and the NASA Shuttle tracking sources produce 
state vectors using the M50 ECI reference frame, and GPS state vectors are provided in 
the WGS 84 Earth-Centered Earth-Fixed (ECEF) reference frame. As shown in Figure 
2.1, the M50 frame is fixed in inertial space while the WGS 84 frame rotates with the 
Earth. In order to compare state vectors from both systems it is necessary to perform a 
transformation between the two reference frames. This transformation requires an 
intermediate step of converting the data into the J2000 ECI reference frame. For this 
reason comparison of state vectors from the two sources was conducted in the J2000 
frame. 




Figure 2.1. ECI and ECEF Coordinate Systems. From Larson and Wertz, 1992, p. 95. 



4. Up-East-North Reference Frame 

GPS velocity data is provided in the Up-East-North (UEN) reference frame which 
has its origin fixed with the spacecraft as shown in Figure 2.2. The UEN system can be 
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visualized using the idea of the Shuttle flying with its belly facing the Earth much like the 
attitude an aircraft would assume while flying above the surface of the Earth. The 
Shuttle’s nose would point along the velocity vector. Up would be equivalent to altitude 
above the ground. East would point straight ahead from the nose and North would point 
in the direction of the left wing tip. Technically, the Up axis is in the radial direction from 
the center of the Earth. The East axis is in the local orbital downtrack direction and 
parallels the velocity vector. The North axis is in the crosstrack direction completing a 
right-handed coordinate system and parallels the orbit normal. 



Z 




RADIAL 



Figure 2.2. Up-East -North Reference Frame. 

In order to perform the transformation from the UEN reference frame to the WGS 
84 reference frame two Euler angle rotations must be performed. A rotation about the Y- 
axis by +8 followed by a rotation about the Z-axis by a value of -y is required as shown in 
Equations 2.2, 2.3 and 2.4. 
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(2.4) 



This transformation is easily realizable in Fortran and was applied to GPS state 
vector velocity components. 



5. Time 

International Atomic Time (TAI) compares several processes of physics such as 
cesium frequency standards and hydrogen masers to form a standard timescale that can be 
used to uniquely identify the instance of time at which an event occurs on Earth. 
(Seidelmann, 1992, p. 2) Universal Time (UT) is based upon the Earth’s rotation with 
respect to the Sun. Since the rotation of the Earth is affected by irregular forces, UT is 
irregular with respect to TAI. The difference between TAI and UT is increasing irregularly 
by approximately one second every eighteen months, and this difference is always 
maintained as an integral number of seconds known as leap seconds. Coordinated 
Universal Time (UTC) is an atomic timescale that is kept in close agreement with UT 
through the addition of these leap seconds. (Seidelmann, 1992, pp. 6-7) UTC is the 
timescale used by GPS and is provided in the receiver navigation solution. 

Greenwich Mean Time (GMT) is the basis of civil time for the United Kingdom 
and is subsequently related to UTC; however, in navigation terminology GMT has been 
used as UT. As a result, the term GMT is somewhat ambiguous. (Seidelmann, 1992, p. 7) 
All Shuttle tracking data uses GMT as its timescale. In this specific case, GMT is 
equivalent to UTC. (Nikolaidis, e-mail, June 1996) In order to facilitate the use of the 
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coordinate frame transformation program and STK, all times were converted to Mission 
Elapsed Time (MET) in seconds which starts at the time of launch. 

B. SHUTTLE TRACKING 

NASA uses two primary sources for Shuttle tracking data, ground sites using C- 
band radar and Tracking and Data Relay Satellites (TDRS). Several of the NASA 
tracking sources are listed in Table 2.1. C-band radar sites use range and angle data and 
TDRS uses doppler measurements to create Shuttle state vectors. Normally two to four 
state vectors are generated per orbit, although, more vectors are generated during periods 
of high activity such as rendezvous, translational maneuvers and LMU alignments. With 
the exclusion of rendezvous, the three sigma position accuracy of the state vectors is 200 
m in the radial direction, 450 m in the downtrack direction and 200 m in the crosstrack 
direction. A three sigma velocity accuracy of 0.45 m/s in the radial direction, 0.20 m/s in 
the downtrack direction and 0.25 m/s in the crosstrack direction is typically achieved. 
During quiescent periods, three sigma accuracies of 100 m and 0.30 m/s radially, 250 m 
and 0.10 m/s downtrack and 100m and 0.15 m/s crosstrack may be obtained. During 
rendezvous and other highly active periods errors on the order of several kilometers may 
be experienced. (RSOC, 1996, p 1) 



Source 


Identifier 


East TDRS 


ESTR 


West TDRS 


WSTR 


Ascension Island 


ASCC 


Bermuda 


BDQC 


Kwajalein Island 


KMTC 


Wallops Island 


WLRC 



Table 2.1. NASA Shuttle Tracking Sources. 
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m. DATA FILES, TOOLS AND PROCEDURE 



A. DATA FILES 



L GPS Navigation Solutions 

GPS data files from the 3M receiver for each flight day of STS-69, Julian days 250 
through 260, were made available by the Johnson Space Center (JSC). There was an 
average of 60 navigation solution files per day and a total of 660 files. Each file had from 
800 to 3300 state vectors with fixes obtained generally every second, although every day 
had some period of outage with no GPS data. Individual state vector entries contained: 
GPS system time; Coordinated Universal Time (UTC); position in XYZ WGS 84 ECEF 
coordinates (m); velocity in UEN coordinates (m/s); position in latitude, longitude and 
altitude; pitch, roll and heading information; acceleration in UEN coordinates; and receiver 
health and status information. A typical GPS data file is shown in Figure 3.1. 

2. NAVG-11 State Vectors 

Rockwell Space Operations Company (RSOC) provided a NAVG-11 Shuttle 
Navigation Postflight Product which included ground navigation history vectors. The 
NAVG-11 product contains the orbiter solution vectors which are processed real-time 
using C-band range and angle data and TDRS doppler data. Each entry contained the 
batch ID, orbit number, batch number, UTC labeled as Greenwich Mean Time (GMT), 
and position and velocity in XYZ M50 ECI coordinates (ft and ft/s). There were 
approximately 430 state vectors for the entire mission in this file. Two to four vectors 
were provided per orbit with the bulk of the data coming from TDRS. (RSOC, 1996, p. 
1) A portion of the NAVG-1 1 file is shown in Figure 3.2. 
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3. 



Propagated State Vectors 



RSOC also provided a second more densely propagated data file with over 1 6,000 
state vectors for the entire mission. The first column of the propagated file is Mission 
Elapsed Time (MET) in hours from launch which was given to be 1995:250:15:09:00.0 
GMT. The program that produces the denser ephemeris propagates from vector to vector 
using a Nystrom-Lear fourth order integrator for both position and velocity. Four 
derivative evaluations are required for each 60 second integration step. The program 
takes into account modeled translational maneuvers. Sun and Moon induced gravitational 
perturbations, drag computed using a Jacchia 1970 atmosphere model, and a 7x7 
gravitational potential model. Figure 3.3 shows a small portion of the propagated data 
file. (Nikolaidis, e-mail, March 1996) 

B. TOOLS 

1. File Transfer Protocol 

All analysis was performed on a Silicon Graphics Indigo II workstation. File 
Transfer Protocol (FTP) provided the only manageable means for retrieving the GPS data 
files from the JSC server. The large number and size of the files made the option of 
copying them to some form of media for transfer unattractive. 

2. Unix Scripts 

Unix scripts, sequences of Unix commands that can be stored in files, were 
employed in editing all the data files and were crucial in the editing of the GPS data files. 
Hand editing all the data files was not feasible, and the Unix scripts proved to be the only 
way to deal with the large quantity of data files by automating the procedure of 
uncompressing and compressing files, grouping files by day, editing files and executing 
Fortran programs. More specifically, the Stream Editor (SED) command made the 
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analysis of the enormous amount of data possible. It provided an automated means for 
removing text headers and troublesome characters such as colons and slashes. 
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3. 



Fortran Programs 



Fortran programs proved to be very versatile and effective in dealing with the 
varying data files. They were used to select the desired state vector elements from the 
input file line entries, convert the input file times to MET in seconds, and write the new 
data files. In the case of the GPS data files, Fortran programs were also used to thin the 
dense GPS data files by a factor of ten and transform GPS UEN velocity data into WGS 
84 coordinates. Fortran routines also converted all the NAVG-11 and propagated data 
from English units to metric units. 

4. Satellite Tool Kit 

Satellite Tool Kit (STK) is an orbit analysis software package developed by 
Analytical Graphics that employs an object paradigm. A Graphical User Interface (GUI) 
is used to manipulate objects such as scenarios, vehicles, facilities, targets, and sensors. 
This thesis required work with only scenarios and vehicles. A vehicle is defined as a 
movable land, sea, air or space object. A scenario is a set of objects that is to be 
visualized. The objects in the scenario are related by time to a map background which 
provides a geographic reference frame. (Analytical Graphics, 1996, pp. 1-2) STK Vehicle 
files for the Shuttle were created using GPS ephemeris data and the propagated file 
ephemeris data and incorporated into scenarios for visualization. 

5. ConvertTool 

The STK math utility convertTool provided a means for transforming data files 
from one reference frame to another. In particular it permitted files to be transformed 
from M50 to J2000, J2000 to M50, ECEF to J2000, and J2000 to ECEF. The program 
accounts for general precession according to IAU 1976 and general nutation according to 
IAU 1980. It does not account for the Earth’s irregular rotation rate or pole wander. 
(Woodbum, e-mail, 1996) 
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6 . 



Matlab 



Matlab is an interactive numerical computation, data analyis and plotting software 
package. The basic Matlab data element is a matrix, and it uses matrix manipulation to 
perform many numeric calculations in a fraction of the time it would take to write and run 
a program in a language such as Fortran which would perform the same function. Matlab 
allows the execution of sequences of commands that can be stored in files called M-files. 
Matlab was used extensively to compare GPS data and Shuttle reference track data. 

C. PROCEDURE 

Upon undertaking a project with such a large amount of data, procedures needed 
to be established to facilitate the transfer and manipulation of the vast data sets. From the 
beginning, using File Transfer Protocol to retrieve the data, to the end, using Matlab to 
perform major coordinate transformations, the emphasis was on using computer programs 
and tools to simplify all aspects of working with the data. The procedure and tools 
utilized in working with the data is outlined in the flow chart in Figure 3.4. 

1. Obtain Data Files 

File Transfer Protocol (FTP) was employed to retrieve the 660 GPS data files from 
the JSC server. The more manageable NAVG-1 1 and propagated files were obtained on 
floppy disk and 8 mm tape, respectively. Grouping the files by flight day provided the 
most efficient way of handling the data. The best results were achieved by plotting each 
day’s worth of data in both STK and Matlab. 

2. Edit GPS Files with Unix Scripts and Fortran Program 
a. Automation with Unix Scripts 

The next step was to edit and manipulate the files using Unix scripts and 
Fortran programs. The following Unix script, “thinfiles,” was invoked to execute the 
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Figure 3.4. Procedure Flow Chart. 
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UEN Transformation Run STK Scenarios 

•uen251.m ^ ‘Equidistant cylindrical plot 

•uen252.m ‘Perspective view plot 





scripts and the Fortran program that would manipulate all 660 GPS data files for flight 
days 250 through 260. 
thinfiles 

fileOscript 
file 1 script 
file2script 
file3 script 
file4 script 
file5 script 
file6script 
file7script 
file8 script 
file9script 
file60script 

Each file#script in “thinfiles” executed the “edit 1 script” script for all the GPS data files in 
each flight day where # was the last digit in the Julian date. 
file 1 script 

edit 1 script 251001 
edit 1 script 251002 
edit 1 script 251003 
edit 1 script 251004 



edit 1 script 251086 

The “edit 1 script” Unix script was the most important script. It uncompressed and 
compressed the raw data files, ran the stream editor command, executed the Fortran 
program, invoked convertTool and created the new edited data files. 



edit 1 script 



uncompress Sl.sol 

sed -f sedl script Sl.sol > dat.edt 

a. out 

cp out.ecf Sl.ecf 
convertTool <out.ecf> $l.j20 
rm dat.edt 
rm out.ecf 
compress Sl.sol 



- uncompresses raw data files 

- runs stream editor and creates data file 

- executes Fortran program 

- copies Fortran program output to new file 

- invokes convertTool program 

- removes intermediary data file 

- removes intermediary data file 

- compresses raw data file for storage 
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The stream editor command eliminated countless hours of hand editing for the raw 
GPS data files by providing an automated way of removing unwanted file header lines and 
troublesome colons and slashes. The script “sedl script” is called by the SED command 
and contains the editing commands that were desired to be performed. 

sedl script 

/GPS/d - deletes file header line which contains the character string GPS 

s/:/ /g - replaces all colons with spaces 

s ! /! ! g - replaces all slashes with spaces 

b. File Editing with Fortran Program 

The Fortran program “thingps” removed the desired state vector elements 

from the edited input file, converted time to MET in seconds, transformed GPS UEN 

velocity data into WGS 84 coordinates, created the new data files, wrote file header 

information required to execute convertTool and thinned the GPS data files by a factor of 

ten. An example of the output from “thingps” is shown in Figure 3.5. 

thingps.f 
program thingps 

implicit none - declare variables 

double precision x,y,z,velx,vely,velz,velup,veleast,velnorth 
double precision gps 
integer yr,doy,hr,min,sec,mo,day 

integer year,dayoyear,month,dayomon,hour,minute, second 
integer tmission,secs 
integer i,count,il 

open(5,file='dat.edt',status='old') - open edited data file 

rewind(5) 

open(6,file- out. ecf, status- unknown') - open new data file 

rewind(6) 

year= 1995 - data required for convertTool 

dayoyear=250 header and MET conversion 

month=09 

dayomon=07 

hour=15 

minute=9 

second=0 

tmission=dayoyear*24*3600+hour*3600+minute*60+second 
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write(*, 1)3 - convertTool header information 

1 format(il) 

write(*,2) year,month,dayomon,hour,minute,second 

2 format(i4,5(5x,i2)) 

count=0 

do 100 i= 1,90000 

read (5, *,end=200) gps,yr,doy,hr,min,sec,x,y,z, - read desired data 

% veleast,velnorth,velup 

secs=(doy * 24*3 600+hr *3 600+min* 60+sec)-tmission -MET conversion 

call uen(x,y,z,velup,veleast,velnorth,velx,vely,velz) 

il = i/10* 10 - thin data files by factor of ten 

if (i.eq.il) then 

write(6,3 ) secs,x,y,z,velx, vely,velz 

end if 

count=count+l 
100 continue 

200 continue 

close(5) 
close(6) 

3 format(i7,6(2x,fl6.6» 

end 

subroutine uen(x,y,z,velup,veleast,velnorth,velx,vely,velz) - velocity conversion 
implicit none 

double precision x,y,z,velup,veleast,velnorth 
double precision R,gamma,d,delta,velx,vely,velz 

R=(x**2+y**2+z**2)**0.5 - elements required for conversion 

gamma=atan2(y,x) (see UEN reference frame section) 

d=(x**2+Y**2)**0.5 
delta=atan2(z,d) 

velx=cos(gamma)*cos(delta)*velup-sin(gamma)*veleast-cos(gamma) -velocity 
% *sin(delta)*velnorth components 

vely=sin(gamma)*cos(delta) :| 'velup+cos(ganima)*veleast- 
% sin(gamma)*sin(delta)*velnorth 

velz=sin(delta) * velup+cos(delta) * velnorth 



return 

end 
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a Reference Frame Transformation Using ConvertTool 

The convertTool program transformed the Fortran output file “out.ecf’ 



from ECEF coordinates to J2000 coordinates producing a file with the extension “j20”. 
The “out.ecf’ header “3” corresponds to this conversion to J2000, and a “1” would 
correspond to a conversion from M50 coordinates to J2000 coordinates. ConvertTool 
also requires a start time, or epoch, in the file header when converting to or from ECEF 
coordinates. This start time was selected to be the time of launch in order to comply with 
the MET convention. Figure 3.6 shows an example output file from convertTool in J2000 
coordinates. 

3. Group GPS Files by Day 

After editing the raw GPS data files and creating the new files with time in MET 
and position and velocity in J2000 XYZ coordinates, the Unix script “daysscript” was 
utilized to group the data files by day. First, empty files were created for each day and 
labeled “gpsday#,” where # was equal to the last digit of the Julian date. Next, 
“daysscript” moved the empty files to a temporary file, “gps.j20,” for editing by the 
“append#script.” The “Append#scipt” invoked the script “dataappend” to append 
sequentially all the files for a flight day to create a single file. Finally, the edited files were 
moved back to the appropriate “gpsday#” file. 



Daysscript 



append 1 script 



mv gps.j20 gpsday0.j20 
mv gpsdayl ,j20 gps.j20 
append 1 script 



mv gpsday0.j20 gps.j20 
appendOscript 



dataappend 251001 
dataappend 251002 
dataappend 251003 



mv gps.j20 gpsdayl .j20 



dataappend 251086 



mv gpsday60.j20 gps.j20 

append60script 

mv gps.j20 gpsday60.j20 



dataappend 
cp $l.j20 data.j20 
cat data.j20 » gps.j20 
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4. 



Edit Propagated Data File 



The propagated data file required much less manipulation. A Fortran program, 
“edtprop,” was used to add the header for convertTool, convert the input file times to 
MET in seconds and convert the data from English units to metric units. A copy of the 
program is in Appendix A. ConvertTool was then used to transform the file from M50 
coordinates to J2000 coordinates. Finally, the file of over 1 6,000 state vectors was hand 
edited to break it up by flight day. 

5. Create STK Vehicle Files 

STK Vehicle files for each flight day were created using the Shuttle STK vehicle 
library file as a guide. A vehicle file was created for each day from each data source. GPS 
ephemeris data and the propagated file ephemeris data were each inserted into a model 
Shuttle vehicle file in ECI J2000 coordinates. A scenario was created for every flight day 
which included a vehicle generated from GPS ephemeris data and one generated from the 
propagated ephemeris data. An example vehicle file is shown in Appendix C. The 
scenarios were animated so that both vehicles could be viewed simultaneously to better 
visualize the similarities and differences between the tracks. 

6. Matlab Comparison 

Matlab was employed for a detailed comparison of the GPS and propagated data 
files for each day. The Matlab M-file, “diffday#.m” was executed to plot J2000 XYZ 
position and velocity vector components and the position and velocity vector magnitudes 
for both data sets. A plot of the data points in three dimensions was also obtained using 
“orbplt#.m” to assist in visualizing the orbit that was generated by the GPS data. 
Appendices D and E contain examples of these M-files. 

Since the propagated data and GPS data could not be compared point for point, it 
was difficult to obtain hard numbers for comparison of GPS and reference track data. In 
an effort to match the data sets point for point, the process of editing the GPS data files 
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for days 25 1 and 252 was repeated with the exception of thinning the data by a factor of 
ten. The NAVG-11 file was hand edited to remove headers and unnecessary columns, 
and then run through the stream editor to remove colons. The Fortran program “edtsv” 
shown in Appendix B was then executed to convert state vector times to MET and 
English units to metric units. ConvertTool was invoked to transform the file from M50 
coordinates to J2000 coordinates. GPS and NAVG-11 state vectors were then hand 
matched for days 251 and 252. A Matlab M-file was used to compare position and 
velocity vector components, position and velocity vector magnitudes, and Root Mean 
Square (RMS) differences between the GPS data and the NAVG-1 1 data. A copy of the 
Matlab script “diffmat#.m” is shown is Appendix F. 
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IV. RESULTS 



GPS data was available for every flight day of STS-69; however, GPS data was 
not available for the entire duration of each flight day. The percentage of each day for 
which GPS data was available for comparison with the propagated data varied greatly 
from day to day. Day 258 was covered most sparsely by GPS with data being available 
only 12 pecent of the day. In contrast, GPS data was available for 93 percent of day 255. 
The average percentage of the day during the mission for which GPS data was available 
was 65 percent. These results are shown in Figure 4.1. 




251 252 253 254 255 256 257 258 259 260 



Flight Day 

Figure 4.1. Percentage of each day for which GPS data was available. 



A. MATLAB COMPARISON RESULTS 



1. Day 251 

Day 251 data was typical of the results achieved during the mission. Figure 4.2 at 
the end of the section displays how dense the GPS data is relative to the propagated file 
reference track. Of note is the fact that the GPS data, displayed as asterisks, has already 



33 



been thinned by a factor of ten. The propagated file is derived from the original state 
vectors which come from several sources as shown in Figure 4.3. The sources displayed 
trends when compared with the GPS data, and these trends appeared to be diurnal in some 
cases. The list of tracking source identifiers is shown in Table 4.1. 



Source 


Identifier 


East TDRS 


ESTR 


West TDRS 


WSTR 


Ascension Island 


ASCC 


Bermuda 


BDQC 


Kwajalein Island 


KMTC 


Wallops Island 


WLRC 



Table 4.1. NASA Shuttle Tracking Sources. 



A plot of the position vector XYZ components in J2000 coordinates for day 25 1 
shown in Figure 4.4 indicates that the GPS data represented by the heavy line follows very 
closely with the reference track represented by the thinner line. A two hour gap in GPS 
data is seen. This is typical of the gaps in GPS data that were observed every day. A plot 
of the velocity vector in XYZ components in J2000 coordinates in Figure 4.5 also shows 
this gap in GPS data and some spurious GPS data points. Both plots display sinusoidal 
patterns for each of the components which are associated with an inclined circular orbit. 

The position vector magnitude is equivalent to the spacecraft’s distance from the 
center of the Earth, or radius of its orbit, as shown in figure 4.6. Plotting of the position 
vector magnitude in Figure 4.7 shows that the radius is fairly constant as would be 
expected for the Shuttle’s circular orbit. This plot of vehicle altitude was very effective in 
indicating when vehicle maneuvers occurred. The velocity vector magnitude plot also 
indicates that the velocity vector was quite constant as was expected for a circular orbit. 
Both the position and velocity vector magnitude plots show that spurious GPS data points 
were present. 
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The three dimensional plot of the GPS data for day 251 in J2000 coordinates 
shown in Figure 4.8 displays how well the GPS data captured the Shuttle’s circular orbit. 
The fact that the orbit is inclined, and that the GPS data is extremely dense is reflected in 
this plot. Erroneous data points are also visible. 

Point for point comparison of the GPS position vector magnitudes with the 
NAVG-1 1 position vector magnitudes in Figure 4.9 shows an altitude difference of ±150 
m. Figure 4.10 shows that the velocity vector magnitude difference between the two 
data sets was generally ± 1 m/s and did not exceed 4 m/s. The RMS position difference 
between the GPS data and NAVG-1 1 data displayed in Figure 4.11 was primarily between 
1400 and 1550 m, and the RMS velocity difference displayed in Figure 4.12 did not 
exceed 9 m/s. 

2. Typical Data 

The position and velocity vector magnitude plots for day 252 shown in Figure 
4.13 display a sinusoidal pattern associated with what appears to be an elliptical orbit. 
Closer examination of the position vector magnitude plot scale shows that the peak to 
peak difference associated with apogee and perigee is approximately 5000 meters. Gaps 
in GPS data and spurious data points are also experienced. Day 255 position and velocity 
vector magnitude plots displayed in Figure 4.14 are similar to those for day 252; however, 
the errant data points are more severe. 

GPS data for day 258 was very sparse with only 12 percent of the flight day 
covered by GPS data. The sparse amount of data allowed a detailed look at both data sets 
as shown in Figure 4.15 and displayed the imperfect nature of an actual vehicle orbiting 
the Earth. Although there were a few spurious data points, the GPS data followed the 
reference track so closely that it was difficult to discern the GPS data from the reference 
track. 
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3. 



Maneuvers 



The position vector magnitude plot proved extremely useful in detecting Shuttle 
maneuvers. The day 253 position vector magnitude plot shown in Figure 4.16 displays a 
gap in GPS data just prior to the maneuver. A small gap in GPS data is experienced at the 
beginning of the maneuver which involves an altitude change of 10,000 m; however, the 
GPS track follows the reference track through the rest of the maneuver. Figure 4. 17 also 
displays an example of GPS tracking through a maneuver. This maneuver is more severe 
involving a change in altitude of 60,000 m. A loss of GPS data is experienced in the 
middle of the maneuver with GPS data being reacquired before completion of the 
maneuver. 

4. Errant GPS Data 

There are several examples of errant GPS data over the course of the eleven flight 
days of STS-69. One of the best examples occurs on day 250 and is shown in Figure 4.18 
and Figure 4.19. Figure 4.18 is a three dimensional plot of the GPS data for day 250. 
This plot shows a segment of GPS data that fails to conform to the Shuttle’s circular orbit. 
This period of bad data can be clearly identified in the position vector component plot 
shown in Figure 4.19. At the start of the day’s data the GPS data does not follow the 
reference track. After a period of approximately 1300 seconds, the GPS data changes 
drastically and matches the reference track. 

GPS data for day 259 has an interesting anomaly that is exhibited in the position 
vector magnitude plot in Figure 4.20. A segment of GPS data appears to be disjointed 
from the rest of the data. In the three dimensional plot of Figure 4.21 this period of errant 
data appears as a separate orbit at a different inclination. Several spurious data points are 
also present, particularly in the velocity vector magnitude plot. 

The position and velocity vector magnitude plots presented the most dramatic 
results of errant data. The Day 260 position and velocity vector magnitude plot in Figure 
4.22 is an example of this observation. The last 5000 seconds of position and velocity 
data for the Day 260 GPS data varies significantly from the reference track. This variance 
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in the data includes as much as a 20 km difference in altitude and a 1.7 km/s difference in 
velocity. 
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Figure 4.2. Density of GPS Data. 
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Figure 4.4. Day 25 1 Position Vector Components in J2000 Coordinates. 
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Figure 4.5. 



Day 251 Velocity Vector Components in J2000 Coordinates. 
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Figure 4.7. Day 251 Position and Velocity Vector Magnitude. 
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Figure 4.9. Day 251 Position Vector Magnitude Difference. 
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Figure 4.10. Day 251 Velocity Vector Magnitude Difference. 
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Figure 4.11. Day 251 Position RMS Difference. 
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Figure 4.13. Day 252 Position and Velocity Vector Magnitude. 
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Figure 4.14. Day 255 Position and Velocity Vector Magnitude. 
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Figure 4.15. Day 258 Position and Velocity Vector Magnitude. 
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Figure 4.16. Day 253 Position and Velocity Vector Magnitude. 
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Figure 4.17. Day 257 Position and Velocity Vector Magnitude. 
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Figure 4.19. Day 250 Position and Velocity Vector Magnitude. 
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Figure 4.20, Day 259 Position and Velocity Vector Magnitude. 
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Figure 4.22. Day 260 Position and Velocity Vector Magnitude. 
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B. STK RESULTS 



STK proved to be a valuable tool for visualizing the differences between the GPS 
vehicle tracks and the propagated vehicle tracks. The GPS vehicle is labeled 
“STS69gpsday#,” and the propagated file vehicle is labeled “STS69propday#” where # 
corresponds to the last digit of the Julian date. The circle surrounding the vehicle 
represents the vehicle’s field of view as seen by a horizon sensor mounted on the Shuttle. 
This field of view circle was helpful in discerning changes in altitude because it expanded 
as the vehicle’s altitude increased. 

The ground tracks displayed in Figure 4.23 show an example of the extreme 
differences that were occasionally encountered between the two data sets. This example 
occurs on day 250 and provides another look at the discrepancy observed in the Matlab 
plots for this day shown in Figure 4.18 and Figure 4.19. The two tracks eventually did 
coincide as shown in Figure 4.24. When the vehicle tracks coincided in this manner it was 
difficult to distinguish between them because they were plotted directly on top of each 
other. 

Figure 4.25 is more typical of the STK plots obtained. This plot shows the entire 
ground track for day 25 1 . The GPS and propagated ground tracks overlay directly on top 
of each other, and the differences between the two vehicles are indistinguishable because 
they are also plotted on top of each other. Figure 4.26 displays a variance in the data for 
day 251 showing the divergent tracks of the two vehicles. Each day’s data analysis 
showed at least one anomaly similar to this one. The perspective shot in Figure 4.27 
provides a better view of the difference between the two tracks. Figure 4.28 shows a 
similar error for day 253. The expanded field of view for the GPS vehicle indicates that 
there is a difference in altitude between the two vehicles. 

A perspective view is shown in Figure 4.29 and displays how well the vehicles and 
orbits for day 255 match each other. Figure 4.30 is a close-up of the same picture and 
exhibits how closely the vehicles and two sets of orbits overlay on top of each other. 
Figure 4.31 further zooms in and reveals that the vehicles are separated, with the GPS 
track trailing the propagated track by a small distance. 
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Several possibilities exist as to why the GPS data and reference track data do not 
match directly. One possibility is that the error is caused by the transformation between 
reference frames that does not take into account pole wander or the irregular rotation rate 
of the Earth. Another possibility is a spurious data point from the GPS receiver. While 
these factors may contribute to the difference between the data sets, the largest source of 
error is most probably the fact that the files for the GPS track and the files for the 
reference track each had times that were rounded to the nearest whole second. With the 
Shuttle traveling at over 7 km/s this can result in considerable position differences between 
the two data sets particularly in the downtrack direction. 

Figure 4.32 shows a GPS position data point at the center of a box whose sides 
were determined based on the velocity of the Shuttle in the X, Y and Z directions for a 
period of one second at a particular instant on day 25 1 . This error box is approximately 
6000 m long in the X direction, 2500 m wide in the Y direction and 1200 m high in the Z 
direction. As expected, the reference track data point for the same instant fell within the 
box. Figure 4.33 shows a plot of the differences between the GPS data and the NAVG-1 1 
state vectors for day 25 1 . The majority of the differences between the data sets lies within 
the box which encloses the error which may have been induced by a frill second of 
timescale inaccuracy. 

The most extreme example of divergence between the GPS and propagated tracks 
occurred on day 256 and is shown in Figure 4.34. In this plot the perspective view is from 
150,000 km above the Earth and the small disc at the center is the Earth. Figure 4.35 is 
more indicative of the differences between vehicle tracks that were encountered. Typically 
an error in which the GPS track started to lead or lag the propagated vehicle would occur 
once per day. 

The ground track plot for day 258 shown in Figure 4.36 reflects the sparse data for 
this flight day. The perspective view in Figure 4.37 shows how closely the two vehicle 
tracks match. The close-up view in Figure 4.38 also shows that the GPS track coincides 
with the propagated track. Figure 4.39 reveals that the two tracks are right on top of each 
other. The GPS orbit in this plot appears very smooth while the propagated orbit appears 
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to be segmented and not as smooth as the GPS orbit. This is expected since the GPS orbit 
is composed of considerably more data points. 

The plot for day 259 shown in Figure 4.40 indicates some discrepancies between 
the GPS and propagated vehicle ground tracks. It shows two independent vehicles which 
indicates that the vehicle tracks diverge. Figure 4.41 further shows that the GPS data has 
produced a vehicle in a different orbit at a different inclination. This plot matches the 
Matlab three dimensional GPS orbit plot for this day shown in Figure 4.21 . 
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Figure 4.23. Day 250 STK Divergent Shuttle Vehicle Tracks. 
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Figure 4.25. Day 251 STK Coincident Shuttle Vehicle Tracks. 
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Figure 4.27. Day 251 STK Perspective View of Divergent Shuttle Vehicle Tracks. 
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Figure 4.29. Day 255 STK Perspective View of Coincident Shuttle Vehicle Tracks. 
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Figure 4.31. Day 255 STK Perspective View Close-up of Shuttle Vehicle Tracks. 
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Figure 4.32. Day 251 GPS and Reference Positions in J2000 Coordinates. 
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Figure 4.34. Day 256 STK Perspective View of Divergent Shuttle Vehicle Tracks. 
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Figure 4.35. Day 257 STK Perspective View of Divergent Shuttle Vehicle Tracks. 
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Figure 4.36. Day 258 STK Coincident Shuttle Vehicle Tracks. 
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Figure 4.37. Day 258 STK Perspective View of Coincident Shuttle Vehicle Tracks. 
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Figure 4.39. Day 258 STK Perspective View Close-up of Shuttle Vehicle Tracks 
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Figure 4.41. Day 259 STK Perspective View of Divergent Shuttle Vehicle Tracks. 
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C. SHUTTLE UEN COMPARISON RESULTS 

L UEN transformation 

In an effort to better classify the state vector errors encountered as errors in the 
crosstrack direction, in the downtrack direction, or in altitude, the GPS and NAVG-11 
state vectors with matching times were transformed from J2000 coordinates to the 
Shuttle’s UEN frame. In order to perform the transformation from the J2000 reference 
frame to the Shuttle UEN reference frame as shown in Figure 4.42, two Euler angle 
rotations must be performed. A rotation about the Z-axis by a value of y followed by a 
rotation about the Y-axis by -8 are required as shown in Equations 4.1, 4.2 and 4.3. 
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Figure 4.42. Up-East-North Reference Frame. 
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2. UEN results 

This transformation was implemented using Matlab and is available in Appendix G. 
The following results shown in Tables 4.2 and 4.3 and Figures 4.43, 4.44, 4.45, 4.49, 4.50 
and 4.51 at the end of the chapter were obtained from the day 251 and day 252 state 
vector position differences. 





Vertical (m) 


Downtrack (m) 


Crosstrack (m) 


Average 


12.1 


5.5 


0 


Median 


13.2 


85.7 


0 


Standard Deviation 


53.8 


1461.0 


0 



Table 4.2. Day 251 State Vector Position differences. 





Vertical (m) 


Downtrack (m) 


Crosstrack (m) 


Average 


2.4 


-32.7 


0 


Median 


8.0 


-140.0 


0 


Standard Deviation 


53.8 


1348.9 


0 



Table 4.3. Day 252 State Vector Position differences. 
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The vertical position differences shown in Figure 4.43 and 4.49 were less than the 
expected one sigma value of 75 m for SPS. The values obtained for vertical position 
difference matched the values obtained from the position vector magnitude (altitude) 
difference plot when the data sets were analyzed in the J2000 reference frame. These 
values are within the expected SPS 1 56 m vertical positioning accuracy. 

Although the average downtrack position difference was fairly small for both days, 
the extremely large standard deviation indicates that significant downtrack errors were 
encountered. When traveling at 7 km/s a downtrack error of 1400 m corresponds to a 0.2 
s error in time. It was expected that a constant position error in the downtrack direction 
would be observed. The downtrack position differences displayed in Figure 4.44 and 
Figure 4.50 varied approximately ±1400 m, and no constant offset was discovered. These 
results far exceeded the SPS one sigma horizontal positioning error of 50 m. 

Crosstrack errors were all less than lxl 0' 8 m for both days as shown in Figure 4.45 
and Figure 4.51. This was beyond the precision of the data, and makes the differences 
essentially zero. These results were unexpected and did not match the one sigma 
horizontal position error of 50 m for SPS. Due to unidentified sources of error, no firm 
conclusions can be drawn; however, the downtrack and crosstrack errors seem to indicate 
that most of the position error is in the downtrack direction. 





Vertical (m/s) 


Downtrack (m/s) 


Crosstrack (m/s) 


Average 


0.29 


0.83 


-0.16 


Median 


0.70 


0.11 


-0.13 


Standard Deviation 


4.18 


2.15 


0.96 



Table 4.4. Day 251 State Vector Velocity differences. 





Vertical (m/s) 


Downtrack (m/s) 


Crosstrack (m/s) 


Average 


0.58 


1.13 


-0.06 


Median 


0.38 


0.39 


-0.15 


Standard Deviation 


4.08 


2.34 


1.53 



Table 4.5. Day 252 State Vector Velocity differences. 
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The velocity results from the comparison of state vectors for day 251 and day 252 
are shown in Tables 4.4 and 4.5 and Figures 4.46, 4.47, 4.48, 4.52, 4.53 and 4.54. Both 
days’ results are very similar. In general the average values for vertical and downtrack 
velocity differences exceeded the expected SPS one sigma velocity accuracy of 0.5 m/s, 
and the median values were close to this value. The crosstrack velocity data is well within 
the expected SPS one sigma velocity accuracy. 
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Figure 4.43. Day 251 Vertical Position Difference. 
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Figure 4.44. Day 25 1 Downtrack Position Difference. 
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Figure 4.45. Day 251 Crosstrack Position Difference. 
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Figure 4.46. Day 251 Vertical Velocity Difference. 
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Figure 4.48. Day 25 1 Crosstrack Velocity Difference. 
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Figure 4.49. Day 252 Vertical Position Difference. 
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Figure 4.50. Day 252 Downtrack Position Difference. 
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Figure 4.51. Day 252 Crosstrack Position Difference. 
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Figure 4.52. Day 252 Vertical Velocity Difference. 
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Figure 4.53. Day 252 Downtrack Velocity Difference. 
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Figure 4.54. Day 252 Crosstrack Velocity Difference. 
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V. 



CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

GPS navigation solutions were available for approximately 65 percent of the STS- 
69 mission, and they generally coincided with the reference track. Based on vehicle track 
visualization using STK and based on conversion from J2000 coordinates to a spacecraft 
local UEN reference frame, differences between the GPS data and the reference data 
appear to have occurred predominantly in the downtrack direction. The fact that the two 
data sets were rounded to the nearest whole second probably impacted the downtrack 
differences between them significantly. 

State vector differences between the GPS navigation solutions obtained using the 
Standard Positioning Service and the reference NAVG-11 state vectors for day 251 and 
252 produced RMS position differences between the data sets of about 1500 m. One 
sigma position accuracy of 54 m in the vertical direction and approximately 1400 m in the 
downtrack direction were experienced. Velocity vector magnitude differences during this 
period were generally ± lm/s, with a RMS velocity difference of less than 9 m/s. One 
sigma velocity accuracies of approximately 4.2 m/s in the vertical direction, 2.3 m/s in the 
downtrack direction and 1.5 m/s in the crosstrack direction were experienced. A firm 
conclusion regarding GPS accuracies could not be drawn because all sources of error 
were not identified. 

GPS position and velocity data was generally very good; however, spurious data 
points were present. The errors included in these data points were sometimes very 
significant, and outages in GPS data of approximately two hours appeared to occur at 
least once per day. Despite these errors, GPS appeared to be effective in producing good 
state vector data even during vehicle maneuvers. 

Based on these results GPS appears to be an excellent navigation source for 
providing Shuttle state vector information. Considering that this data can be obtained 
real-time and still match the reference so closely, real-time GPS state vector data appears 
even more valuable. An important caveat must be applied to this very valuable GPS data: 
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another navigation source such as INS must be present to provide a check against 
spurious data points and periods of outage. 

B. RECOMMENDATIONS 

Analysis of state vector data from a Precise Positioning Service GPS receiver 
should be conducted. PPS state vector data should provide considerably better position 
and velocity data which is required during rendezvous and proximity operations. STS-79 
will provide an opportunity to assess PPS GPS receiver navigation solutions. 

Before integration with the Shuttle navigation suite, experimentation with a simple 
filter for the GPS receiver navigation solutions is required. Erroneous GPS state vector 
data was sometimes quite severe. A process such as a Kalman filter is required to remove 
navigation solutions outside a specified tolerance. 

A standard library of code in a commonly accepted language such as Fortran 
should be adopted for frequently required astronautical functions such as conversion 
between reference frames. In particular, a simple uniform routine for transformation 
between the J2000 ECI frame and the WGS-84 ECEF frame using general precession 
according to IAU-1976 and general nutation according to IAU-1980 without considering 
pole wander or the irregular rotation of the Earth is required. Algorithms which allow real 
time processing of data such as GPS navigation solutions should be emphasized. In 
contrast, current references emphasize the use of techniques which require post-flight 
processing due to their requirement of celestial observations. 

GPS system time should be adopted by all tracking sources. Ambiguity caused by 
UTC being rounded to the nearest whole second by the GPS receiver and tracking sources 
led to significant errors, particularly in the downtrack direction. Since GPS system time 
can be obtained to a very high accuracy by both the Shuttle receiver and a receiver at the 
tracking site, GPS time would provide a valuable timescale addition to all tracking data 
sets. 
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APPENDIX A. FORTRAN PROGRAM FOR EDITING PROPAGATED FILE 
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program edtprop 

implicit none - declare variables 

double precision hours,xft,yft,zft,vxft,vyft,vzft 
double precision secs,x,y,z,vx,vy,vz 
integer i,count 

open(5,file='prop.out',status-old') - open data file 

rewind(5) 

open(6,file='prop.m50',status- unknown') - open new data file 

rewind(6) 

write(6,*) 1 - convertTool header to 

transform to J2000 

count=0 

do 100 i= 1,90000 

read (5,*,end=200) hours,xft,yfl,zft,vxft,vyft,vzft - read data 

secs=hours*3 600.0 - covert MET hours to seconds 

x=xft*0.3048 - convert feet to meters 

y=yft*0.3048 

z=zft*0.3048 

vx=vxft*0.3048 

vy=vyft*0.3048 

vz=vzft*0.3048 

write(6, 1 ) secs,x,y,z,vx,vy,vz - write to new data file 

count=count+l 
100 continue 

200 continue 

close(5) 

close(6) 

1 format(fl6.6,6(2x,fl6.6)) 

end 
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APPENDIX B. FORTRAN PROGRAM FOR EDITING NAVG-11 STATE 

VECTOR FILE 
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program edtsv 



implicit none - declare variables 

double precision station,xft,yft,zft,vxft,vyft,vzft 
double precision x,y,z,vx,vy,vz 
integer sv,doy,hr,min,sec 

integer year,dayoyear,month,dayomon,hour,minute,second 
integer tmission,secs 
integer i, count 

open(5,file='sv.edt',status-old') - open data file 

rewind(5) 

open(6,file='sv.m50',status= , unknown') - open new data file 

rewind(6) 

year=1995 
dayoyear=250 
month=09 
dayomon=07 
hour=15 
minute=9 
second=0 

tmission=dayoyear*24*3600+hour*3600+minute*60+second -MET start 

write(*,l) 1 - convertTool header to 

1 format(il) transform to J2000 

count=0 

do 100 i=l,90000 

read (5,*,end=200) station,sv,doy,hr,min,sec, - read desired data 

% xft,yft,zft,vxft,vyft,vzft 

secs=(doy*24*3600+hr*3600+Tnin*60+sec)-tmission -MET conversion 

x=xft*0.3048 - convert feet to meters 

y=yft*0.3048 

z=zft*0.3048 

vx=vxft*0.3048 

vy=vyft*0.3048 

vz=vzft*0.3048 

write(6,2) secs,x,y,z,vx,vy,vz - write to new datafile 



- data required for MET 
conversion 



count=count+l 



105 



1 00 continue 

200 continue 

close(5) 

close(6) 

2 format(i7,6(2x,fl6.6)) 

end 



APPENDIX C. STK VEHICLE FILE 
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stk.v.2.0 



BEGIN Vehicle 



Name 

BEGIN VehiclePath 
BEGIN J4Perturbation 


STS69gpsdayl 


EphemEpoch 


7 Sep 1995 15:09:00.000000 


StartTime 

StopTime 

SemiMajorAxis 

AltitudeOfApogee 

AltitudeOfPerigee 

RadiusOfApogee 

RadiusOfPerigee 

PeriodOfOrbit 

MeanMotion 

Inclination 

Eccentricity 

ArgOfPerigee 

LongAscenNode 

Right Ascension 

True Anomaly 

MeanAnomaly 

TimePastAscenNode 

TimePastPerigee 

TimeStep 

OrbitGranularity 

NumberOfPasses 


7 Sep 1995 15:09:00.00 
19 Sep 1995 00:00:00.00 

6748537.000000 

370400.000000 

370400.000000 

6748537.000000 

6748537.000000 
5517.284751 
0.001139 
28.500000 
0.000000000000 
0.000000 
-151.000000 
249.528639 
0.000000 
0.000000 
0.000000 
0.000000 

60.000000 
0.000000 
0 


SizeShapeType 

AscenNodeType 

AnomalyType 

EllipseType 

END J4Perturbation 

BEGIN PassDefn 


AppPeriAlt 
Longitude 
True Anomaly 
Osculating 


Break 

Latitude 

DisplayFlag 


Ascending 

0.000000 

Both 
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FirstPass 



RangeFirstPass 

RangeLastPass 

DisplayScheme 

ScenarioEpoch 



0 

2147483647 

AllEphemeris 

1 Nov 1992 00:00:00.0 



Passes 
END Passes 
END PassDefti 

END VehiclePath 

BEGIN Ephemeris 

NumberOfEphemerisPoints 17342 

ScenarioEpoch 1 Nov 1992 00:00:00.0 

EphemerisEcfTimePosVel 

89944473.000000 -6290411.999970-10669.639545 2444894.749963 -944.961651 - 
7710.190409 -2382.491604 

89944485.000000 -6300850.499989 -92193.726556 2416095.999986 -847.046094 - 
7710.033323 -2420.407165 

89944522.000000 -6325375.499997 -33 1401.656250 2329124.749954 -591.753766 - 
7704.550020 -2514.460050 



90020994.000000 -5868722.494942 883252.418377 3213828.661367 -1541.348168 - 

7851.489098 -89.026636 

90020994.000000 -5868722.494942 883252.418377 3213828.661367 -1541.348168 - 

7851.489098 -89.026636 

90020994.000000 -5868722.494942 883252.418377 3213828.661367 -1541.348168 - 

7851.489098 -89.026636 

END Ephemeris 

BEGIN Attitude 

BEGIN ECFWLH 

AZIMUTH 0.000000 

END ECFWLH 

END Attitude 
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BEGIN Swath 



ElevationAngle 


0.000000 


HalfAngle 


0.000000 


RepType 


NoSwath 


END Swath 
BEGIN Constraints 


ConstraintMask 


0 


MinRange 


0.000000 


MaxRange 


1000000000.000000 


MinAzimuthAngle 


-90.000000 


MaxAzimuthAngle 


180.000000 


MinElevationAngle 


-45.000000 


MaxElevationAngle 


45.000000 


MinGrazingAngle 


0.000000 


MaxGrazingAngle 


45.000000 


MinGrazingAltitude 


185200.000000 


MaxGrazingAltitude 


1852000.000000 


MinGndElevationAngle 


0.000000 


MaxGndElevationAngle 


: 90.000000 


MinSunElevationAngle 


-90.000000 


MaxSunElevationAngle 


90.000000 


LightingCondition 


5 


END Constraints 
BEGIN Attributes 


MarkerColor 


6 


GroundTrackColor 


6 


SwathColor 


6 


LineStyle 


3 


MarkerStyle 


4 


END Attributes 
BEGIN Graphics 



Inherit 


On 


ShowLabel 


On 


ShowGndTrack 


On 
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END Graphics 
BEGIN Extensions 



BEGIN Desc 
END Desc 

BEGIN Group 
END Group 

BEGIN Eclipse 

Penumbra Off 

PenumbraColor 7 

PenumbraLineStyle 0 
PenumbraLineWidth 1 

Sunlight Off 

SunlightColor 7 

SunlightLineStyle 0 

SunlightLineWidth 1 

Umbra Off 

UmbraColor 7 

UmbraLineStyle 0 

UmbraLineWidth 1 

END Eclipse 

BEGIN Exclusion 
END Exclusion 

END Extensions 

BEGIN SubObjects 

Class Sensor 

Horizon 

END Class 
END SubObjects 



END Vehicle 



APPENDIX D. MATLAB M-FELE FOR COMPARING GPS AND 

PROPAGATED DATA 
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% Program: diffdayl.m 
% Jim Jones 
% 25 May 96 

% This program displays x,y and z components of GPS and reference 
% track position and velocity data for one day of data. It also 
% displays the magnitude of the position vector, or altitude, 

% from the center of the earth and the magnitude of the velocity 
% vector. 

% Load the GPS data file for one day 
load gpsdayl.j20 -ascii 
gpsin=gpsdayl; 

% Create time, position component and velocity component vectors. 

gpssec=gpsin(:,l); 

gpsx=gpsin(:,2); 

gpsy=gpsin(:,3); 

gpsz=gpsin(:,4); 

gpsxdot=gpsin(:,5); 

gpsydot=gpsin(: ,6); 

gpszdot=gpsin(:,7); 

% Load the reference track data file for one day 

load propdayl .j20 -ascii 

propin=propdayl; 

% Create time, position component and velocity component vectors. 

propsec=propin(: , 1 ); 
propx=propin(:,2); 
propy=propin(:,3); 
propz=propin(:,4); 
propxdot=propin(: ,5); 
propydot=propin( : ,6); 
propzdot^ropi^: ,7); 
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% Plot position components vs. time 

figure(l) 

subplot(3,l,l) 

plot(gpssec,gpsx,'w*',propsec,propx,'w') 
xlabel('Day 251 Mission Seconds') 
ylabel('X Position (m)') 
grid 



subplot(3,l,2) 

plot(gpssec,gpsy,'w*',propsec,propy,'w') 
xlabel(T)ay 251 Mission Seconds') 
ylabelfY Position (m)') 
grid 

subplot(3,l,3) 

plot(gpssec,gpsz,'w*',propsec,propz,'w') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Z Position (m)') 
tol=-l; 

legend('GPS','Reference',tol) 

grid 

orient tall 

%print_orig -dgifB 251pos.gif 
print 

% Plot velocity components vs. time 

figure(2) 

subplot(3,l,l) 

plot(gpssec,gpsxdot,'w*',propsec,propxdot,'w') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('X Velocity (m/s)') 
grid 

subplot(3,l,2) 

plot(gpssec,gpsydot,’w*',propsec,propydot,'w’) 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Y Velocity (m/s)') 
grid 
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subplot(3,l,3) 

plot(gpssec,gpszdot,'w*',propsec,propzdot,'w') 
xlabel(T)ay 25 1 Mission Seconds') 
ylabel('Z Velocity (m/s)’) 
tol=-l; 

legend('GPS7Reference',tol) 

grid 

orient tall 

%print_orig -dgifS 251vel.gif 
print 

% Calculate position vector magnitudes (altitude) 

gpsr=sqrt(gpsx. A 2 + gpsy. A 2 + gpsz. A 2); 
propr=sqrt(propx A 2 + propy. A 2 + propz. A 2); 

% Calculate velocity vector magnitude 

gpsv=sqrt(gpsxdot A 2 + gpsydot A 2 + gpszdot. A 2); 
propv=sqrt(propxdot A 2 + propydot A 2 + propzdot. A 2); 



% Plot position and velocity vector magnitudes 

figure(3) 

subplot(2,l,l) 

plot(gpssec,gpsr,'w*',propsec,propr,'w') 

xlabel(T)ay 25 1 Mission Seconds') 

ylabel('Position Vector Magnitude (Vehicle Altitude) (m)') 

tol=-l 

legend('GP S'j'Reference'jtol) 
grid 

subplot(2,l,2) 

plot(gpssec,gpsv,'w*',propsec,propv,'w') 
xlabel(T)ay 25 1 Mission Seconds') 
ylabel('Velocity Vector Magnitude (m/s)') 
tol=-l 

legend('GPS','Reference',tol) 

grid 

orient tall 

%print_orig -dgifS 251mag.gif 
print 

% Plot close up of position data 
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figure(4) 

plot(gpssec,gpsx,'w*',propsec,propx,'wo') 
xlabel(’Day 251 Mission Seconds') 
ylabel('X Position (m)') 
tol=-l 

legend('GPS','Reference',tol) 

grid 

axis([(5e4+20) (5e4+220) -0.57e7 -0.47e7]) 

%print_orig -dgifS 251close.gif 

print 
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APPENDIX E. MATLAB M-FILE FOR PLOTTING GPS ORBIT IN 3-D 
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% Program: orbpltl.m 
% Jim Jones 
% 25 May 96 

% This program creates a 3D plot in J2000 coordinates of GPS navigation 
% solutions for one day of flight data. 

% Load GPS ephemeris file in X,Y,Z J2000 coordinates 

load gpsday0.j20 -ascii 

out = gpsdayO; 

% Plot ephemeris from data file 

xl = [0 le7]; 
x2 = [0 0]; 
x3 = [0 0]; 

yl = [0 0]; 
y2 = [0 le7]; 

y3 « [0 0]; 

zl = [0 0]; 
z2 = [0 0]; 
z3 = [0 le7]; 

figure(l) 

plot3(x 1 ,x2,x3 ,'w-',y 1 ,y2,y3 ,'w-',z 1 ,z2,z3 ,'w-',out(: ,2),out(: ,3 ),out(: ,4),'w. ') 

titleCGPS Orbit for Day 250 in J2000 Coordinates') 

xlabel('X Axis (m)') 

ylabel('Y Axis (m)') 

zlabel('Z Axis (m)') 

%print_orig -dgifB 250orb.gif 
print 
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APPENDIX F. MATLAB M-FILE FOR COMPARING GPS AND NAVG-11 
STATE VECTORS WITH MATCHING TIMES 
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% Program: diffmatl.m 
% Jim Jones 
% 25 May 96 

% This program plots differences between GPS navigation solution 
% data and Rockwell state vector data which share the same time. 
% Both data sets are in J2000 coordinates. 



% Load GPS data file 
load gpsmatchl.j20 -ascii 
gpsin=gpsmatchl; 

% Create time, position component and velocity component vectors. 

gpssec=gpsin(:,l); 

gpsx=gpsin(:,2); 

gpsy=gpsin(:,3); 

gpsz=gpsin(:,4); 

gpsxdot=gpsin(:,5); 

gpsydot=gpsin( : ,6); 

gpszdot=gpsin(:,7); 

% Load Rockwell state vector data 

load svmatchl.j20 -ascii 

propin=svmatchl; 

% Create time, position component and velocity component vectors. 

propsec=propin(: , 1 ); 

propx=propin(:,2); 

propy=propin(:,3); 

propz=propin(:,4); 

propxdot=propin(:,5); 

propydot=propin(: ,6); 

propzdot=propin(: ,7); 
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% Identify sources of Rockwell state vectors 
% (Sources corresponding to indices of matrix propin) 


% 


Legend 


estr=[l 3 5 7 9 1 1 16 21 24 27 30 32 34 36]; 


% * 


wstr=[24 6 8 10 13 17 18 22 28 35 37]; 


% 0 


bdqc=[12 15 202631]; 


% X 


wlrc=[14 19 25]; 


% + 


kmtc=[23 29 33]; 


% . 



% Plot position components vs. time 

figure(l) 

subplot(3,l,l) 

plot(gpssec,gpsx,'w*',propsec,propx,'wo') 
xlabel(Day 25 1 Mission Seconds') 
ylabel('X Position (m)') 
grid 

subplot(3,l,2) 

plot(gpssec,gpsy,'w*',propsec,propy,'wo') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Y Position (m)') 
grid 

subplot(3,l,3) 

plot(gpssec,gpsz,'w*',propsec,propz,W) 
xlabel(Day 25 1 Mission Seconds') 
ylabel('Z Position (m)') 
tol=- 1 ; 

legend('GPS','Reference',tol) 

grid 

orient tall 

%print_orig -dgif8 251posm.gif 
print 

% Plot velocity components vs. time 

figure(2) 

subplot(3,l,l) 

plot(gpssec,gpsxdot,'w*',propsec,propxdot,W) 
xlabel(T)ay 25 1 Mission Seconds') 
ylabel('X Velocity (m/s)') 
grid 
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subplot(3,l,2) 

plot(gpssec,gpsydot,'w*',propsec,propydot,'wo') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Y Velocity (m/s)') 
grid 

subplot(3,l,3) 

plot(gpssec,gpszdot,'w*',propsec,propzdot,'wo') 
xlabel('Day 251 Mission Seconds') 
ylabel('Z Velocity (m/s)') 
toI=-l, 

legend('GPSVReference',tol) 

grid 

orient tall 

%print_orig -dgifB 251velm.gif 
print 

% Calculate position vector magnitudes (altitude) 

gpsr=sqrt(gpsx. A 2 + gpsy. A 2 + gpsz. A 2); 
propn=sqrt(propx. A 2 + propy. A 2 + propz. A 2); 

% Calculate velocity vector magnitude 

gpsv=sqrt(gpsxdot A 2 + gpsydot A 2 + gpszdot A 2); 
propv=sqrt(propxdot. A 2 + propydot. A 2 + propzdot A 2); 



% Plot position vector and velocity magnitudes 

figure(3) 

subplot(2,l,l) 

plot(gpssec,gpsr,'w*',propsec,propr,'wo') 

xlabel(T)ay 25 1 Mission Seconds') 

ylabel('Position Vector Magnitude (Vehicle Altitude) (m)') 

grid 

subplot(2,l,2) 

plot(gpssec,gpsv,'w*',propsec,propv,'wo') 
xlabel(T)ay 25 1 Mission Seconds') 
ylabel('Velocity Vector Magnitude (m/s)') 
tol=-l 

legendCOPS^Tleference^tol) 

grid 

orient tall 

%print_orig -dgifS 251magm.gif 
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print 

% Calculate position component differences 



posdiffx=gpsx-propx; 

posdiffy=gpsy-propy; 

posdiffz=gpsz-propz; 

% Plot position component differences 

% Data source legend 

% estr- * 

% wstr - o 
% bdqc - x 
% wire - + 

% kmte - . 



figure(4) 

plot(propsec(estr),posdiffx(estr),'w*',... 

propsec(wstr),posdiffx(wstr),'wo’,propsec(bdqc),posdiffx(bdqc),'wx',. 
propsec( wire), posdiffx( wire), 'w+',propsec(kmtc),posdiffx(kmtc),'w.') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('X Position Difference (m)') 
grid 
tol=- 1 ; 

legend('ESTR','WSTR', l BDQC','WLRC','KMTC',tol) 

%print_orig -dgifB 251difxm.gif 
print 

figure(5) 

plot(propsec(estr),posdiffy(estr),'w*',. . . 

propsec(wstr),posdiffy(wstr),'wo',propsec(bdqc),posdiffy(bdqc),'wx',. 
propsec(wlrc),posdiffy(wlrc),'w+',propsec(kmtc),posdiffy(kmtc),'w.') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Y Position Difference (m)') 
grid 
tol=-l; 

legend('ESTR',’WSTR','BDQC’,’WLRC','KMTC’,tol) 

%print_orig -dgifB 251difym.gif 
print 
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figure(6) 

plotCpropse^estrXposdiff^estrX'w*',. . . 

propsec(wstr),posdiffz(wstr),'wo',propsec(bdqc),posdiffz(bdqc),'wx', 
propsec(wlrc),posdiffz(wlrc),'w+',propsec(kmtc),posdiffz(kmtc),'w.') 
xlabel('Day 251 Mission Seconds') 
ylabel('Z Position Difference (m)') 
grid 
tol=-l; 

legend( r ESTR', r WSTR', r BDQCVWLRCVKMTC',tol) 

%print_orig -dgifS 251difzm.gif 
print 



% Calculate velocity component differences 

veldiffx=gpsxdot-propxdot; 

veldiffy=gpsydot-propydot; 

veldiffz=gpszdot-propzdot; 

% Plot velocity component differences 

figure(7) 

plot(propsec(estr),veldiffx(estr),'w*',... 

propsec(wstr),veldiffx(wstr),'wo',propsec(bdqc),veldiffx(bdqc),'wx',.. 
propsec(wlrc),veldiffx(wlrc),'w+',propsec(kmtc),veldiffx(kmtc),'w.') 
xlabel(T)ay 25 1 Mission Seconds') 
ylabel('X Velocity Difference (m/s)') 
grid 
tol=-l; 

legend('ESTR', , WSTR','BDQC','WLRC', , KMTC',tol) 

%print_orig -dgifB 251dfxdm.gif 
print 

figure(8) 

plot(propsec(estr),veldiffy(estr),'w*',. . . 

propsec(wstr),veldiffy(wstr),'wo',propsec(bdqc),veldiflfy(bdqc),'wx',.. 
propsec(wlrc),veldiffy(wlrc),'w+',propsec(kmtc),veldiffy(kmtc),'w.') 
xlabel(T)ay 25 1 Mission Seconds') 
ylabelCY Velocity Difference (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','KMTC',tol) 

%print_orig -dgifB 251dfydm.gif 
print 
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figure(9) 

plot(propsec(estr),veldiffz(estr),'w*',... 

propsec(wstr),veldiffz(wstr),'wo',propsec(bdqc),veldiffz(bdqc),'wx',... 
propsec(wlrc),veldiffz(wlrc),'w+',propsec(kmtc),veldiffz(kmtc),'w.') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Z Velocity Difference (m/s)') 
grid 
tol— 1 ; 

legend('ESTR','WSTR', r BDQCVWLRC','KMTC',tol) 

%print_orig -dgifS 251dfzdm.gif 
print 



% Calculate position and velocity magnitude differences 

diffr=gpsr-propr, 

diffv=gpsv-propv; 

% Plot position and velocity magnitude differences 
figure(lO) 

plot(propsec(estr),diffr(estr),'w*',... 

propsec(wstr),diffr(wstr),'wo',propsec(bdqc),difif(bdqc),'wx',... 
propsec(wlrc),diflEf(wlrc),'w+',propsec(kmtc),diffr(kmtc),'w.') 
xlabel('Day 25 1 Mission Seconds') 

ylabel('Position Vector Magnitude (Altitude) Difference (m)') 

grid 

tol=- 1 ; 

legend( r ESTR', r WSTR', r BDQCVWLRC','KMTC',tol) 

%print_orig -dgifS 251difrm.gif 
print 

figure(ll) 

plot(propsec(estr),difiV(estr),'w*',... 

propsec(wstr),difiv(wstr),'wo',propsec(bdqc),diffV(bdqc),'wx', . . . 
propsec(wlrc),diffv(wlrc),'w+',propsec(kmtc),dif3V(kmtc),'w.') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Velocity Vector Magnitude Difference (m/s)') 
grid 
tol=-l; 

legend('ESTR', r WSTR','BDQC','WLRC','KMTC',tol) 

%print_orig -dgifB 251difvm.gif 
print 
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% Calculate position and velocity RMS differences 
rmsi=sqrt(posdiffx. A 2 + posdiffy. A 2 + posdiffz. A 2); 
rmsv=sqrt(veldiffx. A 2 + veldiffy. A 2 + veldiffy. A 2); 

% Plot position and velocity RMS differences 
figure(12) 

plot(propsec(estr),rmsr(estr),'w*' v . . 

propsec(wstr),rmsr(wstr), , wo',propsec(bdqc),rmsr(bdqc),’wxV 
propsec(wlrc),rmsr(wlrc),'w+',propsec(kmtc),rrnsr(kmtc),'w.') 
xlabel('Day 251 Mission Seconds') 
ylabelCPosition RMS Difference (m)') 
grid 
tol=-l; 

legend('ESTR', r WSTR'3DQC','WLRC','KMTC',tol) 

%print_orig -dgifB 251drmsm.gif 
print 

figure(13) 

plot(propsec(estr),rmsv(estr),'w :,: ',. . . 

propsec(wstr),rmsv(wstr),'wo',propsec(bdqc),rmsv(bdqc),'wx', 
propsec(wlrc),rmsv(wlrc),'w+',propsec(kmtc) ? rmsv(kmtc),'w.') 
xlabel(Day 251 Mission Seconds') 
ylabel('Velocity RMS Difference (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','KMTC’,tol) 

%print_orig -dgifB 251dvmsm.gif 
print 
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APPENDIX G. MATLAB M-FILE FOR TRANSFORMING MATCHING GPS 
AND NAVG-11 STATE VECTORS TO UEN FRAME FOR COMPARISON 
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% Program: uen251.m 
% Jim Jones 
% 25 June 96 

% This program plots differences between GPS navigation solution 
% data and NAVG-1 1 state vector data which share the same time. 
% Both data sets are in the Shuttle UEN frame. 

% Load the GPS data file for one day 

load gpsmatchl.j20 -ascii 

gpsin=gpsmatchl; 

% Create time, position component and velocity component vectors. 

gpssec=gpsin(:,l); 

gpsx=gpsin(:,2); 

gpsy=gpsin(:,3); 

gpsz=gpsin(:,4); 

gpsxdot=gpsin(:,5); 

gpsydot=gpsin(: , 6), 

gpszdot=gpsin(:,7); 

% Load the reference track data file for one day 

load svmatchl.j20 -ascii 

svin=svmatchl; 

% Create time, position component and velocity component vectors. 

refsec=svin(:,l); 

refx=svin(:,2); 

refy=svin(:,3); 

refz=svin(:,4); 

refxdot=svin(:,5); 

refydot=svin(:,6); 

refzdot=svin(:,7); 
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% Identify sources of NAVG- 11 state vectors 
% (Sources corresponding to indices of matrix svin) 



% 


Legend 


estr=[ 1 3 5 79 11 16 21 24 27 30 32 34 36];% * 


wstr=[2 4 6 8 10 13 17 18 22 28]; 


% 0 


bdqc=[12 15 20 26 31]; 


% X 


wlrc=[14 19 25]; 


% + 


kmtc=[23 29 33]; 


% . 



for i=l:length(gpsx) 

% Elements required to calculate DCM for GPS data 

gpsR=sqrt(gpsx(i) A 2+gpsy(i) A 2+gpsz(i) A 2); 

gpsgamma=atan2(gpsy(i),gpsx(i)); 

gpsd=sqrt(gpsx(i) A 2+gpsy(i) A 2); 

gpsdelta=atan2(gpsz(i),gpsd); 

% DCM components for GPS data 

gpsdcm 1 1 =cos(gpsdelta) * cos(gpsgamma) ; 
gpsdcm 1 2=cos(gpsdelta) * sin(gpsgamma); 
gpsdcm 1 3=sin(gpsdelta); 

gpsdcm2 1 =sin(gpsdelta); 

gpsdcm22=cos(gpsgamma); 

gpsdcm23=0; 

gpsdcm31=-sin(gpsdelta)*cos(gpsgamma); 
gpsdcm3 2=-sin(gpsdelta) * sin(gpsgamma) ; 
gpsdcm3 3 =cos(gpsdelta) ; 

% DCM rows for GPS data 

gpsdcm l=[gpsd cm 1 1 gpsdcm 12 gpsdcm 13]; 
gpsdcm2=[gpsdcm2 1 gpsdcm22 gpsdcm23]; 
gpsdcm3=[gpsdcm3 1 gpsdcm32 gpsdcm33]; 

% 3x3 DCM for GPS data 

gpsdcm=[gpsdcm 1 ; gpsdcm2; gpsdcm3]; 
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% 



Create GPS position and velocity vectors 



gpspos=[gpsx(i); gpsy(i); gpsz(i)]; 
gpsvel=[gpsxdot(i); gpsydot(i); gpszdot(i)]; 

% Perform transformation from J2000 coordinates to UEN 

gpsuenp=gpsdcm*gpspos; 

gpsuenv=gpsdcm*gpsvel; 

% Elements required to calculate DCM for reference data 

refR=sqrt(refx(i) A 2+refy(i) A 2+refz(i) A 2); 

refgamma=atan2(refy(i),refx(i)); 

refd=sqrt(refx(i) A 2+refy(i) A 2); 

refdelta=atan2(refz(i),refd); 

% DCM components for reference data 

refdcml 1 =cos(refdelta)*cos(refgamma); 
refdcm 1 2=cos(refdelta)*sin(refgamma); 
refdcm 1 3 =sin(refdelta); 

refdcm2 1 =sin(refdelta); 

refdcm22=cos(refgamma); 

refdcm23=0; 

refdcm31=-sin(refdelta)*cos(refgamma); 
refdcm3 2=-sin(refdelta)* sin(refgamma); 
refdcm3 3=cos(refdelta); 

% DCM rows for reference data 

refdcm 1 =[refdcm 1 1 refdcm 12 refdcm 13]; 
refdcm2=[refdcm21 refdcm22 refdcm23]; 
refdcm3=[refdcm3 1 refdcm32 refdcm33]; 

% 3x3 DCM for reference data 

refdcm=[refdcml; refdcm2; refdcm3]; 

% Create reference position and velocity vectors 

refpos=[refx(i); refy(i); refz(i)]; 
refvel=[refxdot(i); refydot(i); refzdot(i)]; 
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% 



Perform transformation from J2000 coordinates to UEN 



refuenp=refdcm * refpos; 
refuenv=refd cm * refvel ; 

% Calculate state vector differences 

difp=gpsuenp-refuenp; 

difv=gpsuenv-refuenv; 

vertp(i)=difp(l); 
dt rkp(i)=di^>(2 ) ; 
xtrkp(i)=difjj(3); 

vertv(i)=difv(l); 

dtrkv(i)=difv(2); 

xtrkv(i)=difV(3); 



end 

% Calculate average differences 

vertpavg=mean(vertp); 

dtrkpavg=mean(dtrkp); 

xtrkpavg=mean(xtrkp); 

vertvavg=mean(vertv); 

dtrkvavg=mean(dtrkv); 

xtrkvavg=mean(xtrkv); 

% Calculate median differences 

vertpmed=median(vertp); 

dtrkpmed=median(dtrkp); 

xtrkpmed=median(xtrkp); 

vertvmed=median(vertv); 

dtrkvmed=median(dtrkv); 

xtrkvmed=median(xtrkv); 
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% Calculate standard deviations of differences 



vert pst d=std( vertp) ; 
dtrkpstd=std(dtrkp); 
xtrkpstd=std(xtrkp); 

vertvstd=std(vertv); 
dtrkvstd=std(dtrkv); 
xt rkvst d=st d(xtrkv); 

% Plot position component differences 

figure(l) 

plot(gpssec(estr),vertp(estr),'w*',gpssec(wstr),vertp(wstr),'wo',... 
gpssec(bdqc),vertp(bdqc),'wx'. > gpssec(wlrc). > vertp(wlrc),'w+'. > . . . 
gpssec(kmtc), vertp(kmtcX'w. ') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Vetical Position Difference (m)') 
grid 
tol=- 1 ; 

legend('ESTRVWSTR','BDQC', r WLRC','KMTC',tol) 

print 

figure(2) 

plot(gpssec(estr),dtrkp(estr),'w*',gpssec(wstr),dtrkp(wstr),'wo',... 
gpssec(bdqc),dtrkp(bdqc),'wx',gpssec(wlrc),dtrkp(wlrc),'w+',... 
gpssec(kmtc),dtrkp(kmtc),'w.') 
xlabel('Day 25 1 Mission Seconds') 
ylabel(T)owntrack Position Difference (m)') 
grid 
tol=- 1 ; 

legend('ESTR', r WSTR','BDQC','WLRC','KMTC',tol) 

print 

figure(3) 

plot(gpssec(estr),xtrkp(estr),'w*',gpssec(wstr),xtrkp(wstr),'wo',... 
gpssec(bdqc),xtrkp(bdqc),'wx',gpssec(wlrc),xtrkp(wlrc),'w+',. . . 
gpssec(kmtc),xtrkp(kmtc),'w.') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Crosstrack Position Difference (m)') 
grid 
tol=-l; 

legend( r ESTRVWSTR','BDQC','WLRC','KMTC',tol) 

print 
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% Plot velocity component differences 
figure(4) 

plot(gpssec(estr),vertv(estr),'w*',gpssec(wstr),vertv(wstr),'wo',... 
gpssec(bdqc),vertv(bdqc),'wx',gpssec(wlrc),vertv(wlrc),’w+',... 
gpssec(kmtc),vertv(kmtc),'w. ') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Vetical Velocity Difference (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQCVWLRC','KMTC',tol) 

print 

figure(5) 

plot(gpssec(estr),dtrkv(estr),'w*',gpssec(wstr),dtrkv(wstr),'wo',... 
gpssec(bdqc),dtrkv(bdqc),'wx',gpssec(wlrc),dtrkv(wlrc),'w+',. . . 
gpssec(kmtc),dtrkv(kmtc),'w.') 
xlabel(Day 25 1 Mission Seconds') 
ylabel('Downtrack Velocity Difference (m/s)') 
grid 
tol=- 1 ; 

legend( r ESTR',’WSTR', r BDQCVWLRC','KMTC',tol) 

print 

figure(6) 

plot(gpssec(estr),xtrkv(estr),'w*',gpssec(wstr),xtrkv(wstr),'wo',. . . 
gpssec(bdqc),xtrkv(bdqc),'wx',gpssec(wlrc),xtrkv(wlrc),'w+',... 
gpssec(kmtc),xtrkv(kmtc),'w.') 
xlabel('Day 25 1 Mission Seconds') 
ylabel('Crosstrack Velocity Difference (m/s)') 
grid 
tol=-l; 

legend('ESTR’,'WSTR','BDQC’,'WLRCVKMTC',tol) 

print 
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APPENDIX H. DAY 250 STK PLOTS 
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APPENDIX J. DAY 251 STK PLOTS 
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APPENDIX K. DAY 251 MATLAB PLOTS 
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